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Abstract. Homotopy continuation provides a numerical tool for computing 
the equivalence of a smooth variety in an intersection product. Intersection 
theory provides a theoretical tool for relating the equivalence of a smooth vari- 
ety in an intersection product to the degrees of the Chern classes of the variety. 
A combination of these tools leads to a numerical method for computing the 
degrees of Chern classes of smooth projective varieties in P". We illustrate the 
approach through several worked examples. 



1. Introduction 

In this paper we describe a numerical method for computing the degrees of the 
Chern classes of a smooth projective variety. The main tools that will be used are 
drawn from intersection theory and numerical algebraic geometry. In particular, 
we will be using a refined Bezout's theorem [8l[9l[3T] together with ideas related to 
numerical polynomial algebra such as homotopy continuation, monodromy, and the 
numerical decomposition of zero sets [22j [25l [26l [28j . In addition, we will need to 
numerically determine the degree of a certain residual zero-scheme [Sj El [HI [20l [21] . 
Implementations of many of the computational tools that will be used can be found 
in freely available software packages (e.g. [2l|5l[iail3l[l2[19l[23|30l[32l[34]). The 
algorithms developed in this paper primarily use the software package Bertini with 
help from the software package Macaulay 2 f2l [l2] ■ 

From the generators of the ideal of a smooth projective variety, Chern numbers 
can be computed through a purely symbolic computation. This naturally leads 
one to ask if a numerical approach is useful or needed. A partial answer is that in 
the numerical approach described in this paper, the algorithms work equally well 
whether the generators are sparse or dense and whether they have rational, algebraic 
or transcendental coefficients. This is not the case with purely symbolic methods. 
Furthermore, a surprising amount of information can be extracted from an ideal 
even in situations where the generators have inaccuracies in their coefficients [T51 
1171 123j . Another important feature, that will likely play a dominant role in the 
future, is that homotopy continuation algorithms parallelize well. This will allow 
such algorithms to take advantage of the trend of more processors on a chip. These 
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features lead us to believe that the introduction of a numerical approach is indeed 
useful and will allow experimentation in settings outside the domain of purely 
symbolic methods. 

Acknowledgments: The authors would like to thank Frank-Olaf Schreyer and 
Charles Wampler for stimulating conversations and observations. 

2. INTERSECTION THEORY AND HOMOTOPY CONTINUATION 

In this section, we describe the main tools that will be used from intersection 
theory and from numerical algebraic geometry. Let Xi , . . . , Xr be a set of hyper- 
surfaces in P*". Let Z he a smooth connected component of the scheme defined 
by Xi n . . . n Xr- We present a well known theorem from intersection theory that 
relates the equivalence of Z in the intersection product Xi ■ . . . ■ Xr to a formula 
involving the total Chern classes of certain bundles related to Xi,. . . ,Xr, Z, and 
P''. With mild assumptions on the scheme determined by the hypersurfaces, this 
theorem can be translated into a numerical condition involving the degrees of the 
Chern classes of Z, the degrees of the Xi and the degree of a certain residual 
zero-scheme. Tools from numerical algebraic geometry can be used to compute the 
degree of the residual zero-scheme. By varying the degrees of the Xi, we obtain 
multiple numerical conditions that can be used to uniquely determine the degrees 
of the Chern classes of Z. 

2.1. Intersection Theory. Intersection theory has been primarily developed from 
the two different viewpoints found in [51 [HI [31] ■ Though each viewpoint seems 
to have its own advantages and disadvantages, Leendert van Gastel showed that 
the two are closely related [29]. While cither version would be adequate for our 
purposes, we follow the approach of Fulton-MacPherson as described in [9]. We 
first need to fix the notation and definitions that will be used in this section. 

Definition 2.1. Let Z and X be projective varieties in P*" = PJJ. with Z C X. 

(1) Let A^,{X) — Ak{X) where Ak{X) denotes the Chow group of k- 
dimensional cycles on X modulo rational equivalence. 

(2) For a £ A^,{X), {a}^ denotes the component of a in Ao{X). Let a £ 
Ak{X), a = I^ili where oi, . . . , a™ G Z and Vi, . . . ,V,n ^ X are 
irreducible subvarieties. The degree of a is defined by deg a — '^iCa deg Vi . 

(3) NxP"^ denotes the normal bundle of X inP'^, Tx denotes the tangent bundle 
of X and Tx\z denotes its restriction to Z. 

(4) Given a vector bundle E of rank n on a smooth variety, Ci {E) denotes the 
J*'' Chern class of E and c{E) = 1-1- Ci{E) -I- • • • -I- Cn{E) denotes the total 
Chern class of E. 

(5) The i^^ Chern class of a smooth variety X is defined as Ci{Tx) and 
deg (ci(Tx)) is called a Chern number of X . 

Let Xi , . . . , Xr be hypersurfaces in P'' and let Z be a smooth connected 
component of Xi n . . . fl X^. In [9] Fulton considers the intersection product 
Xi- . . .-Xr e Ao{Xiri. . .nXr) and defines {Xi ■ . . .-XrY G ^o(^), the equivalence 
oi Z ior Xi ■ . . . ■ Xr, as the part of the intersection product supported on Z. The 
following is a simplified version of Proposition 9.1.1 of [9] specialized to our setting. 
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Proposition 2.2. (Fulton) Let Xi, . . . ,Xr be hypersurfaces in and let Z be a 
smooth connected component of Xi fl . . . fl Xr- Let Ni be the restriction of NxJ?^ 
to Z . Then 

(1) (Xi ....-Xrf = {{IiUAmc{TrAzr'c{Tz)},. 

We introduce notation for the purpose of stating a useful corollary of Proposi- 
tion [221 Let m = degXi, let Am\Z = n, and let co, . . . , c„ be the Chern classes of 
Z. The /c"* elementary symmetric function in ni, . . . , will be denoted by (Tfc. In 
other words, 

(70 = 1, cri=^ni, a2 = '^ninj, 0-3 = ^ n.^UjUk, ... 

i i<j i<j<k 

If we let 

j=o V J / 

then we obtain the following identity as a direct corollary of Proposition 12.21 
Corollary 2.3. deg [Xi ■ . . . ■ X,)^ = Y.l=o degc,. 

Proof. Let H be the hyperplane class in P'' and let H\z be its restriction to Z. 
Because c[Tr^\z) = (1 + H\zY+\ 

c{TrAz)-'=jZ{-iy{"^AH\'z- 
3=0 \ J / 

As cycle classes on P'', iVx^P'^ = Xi — UiH. Hence c{Ni) = 1 + niH\z and we get 
that n[^j^(l + c{Ni)) = J22=o '^kH\\. Thus the formula ([T|) can be written as 

fe=0 j"=0 \ J / 

The statement now follows since degQ — deg(cii/|2 *) and the component of 
{Xi ■ . . . ■ X,.)^ in Ao{Z) is given by the sum of terms where k = n — i — j. □ 

The classical version of Bezout's theorem (Proposition 8.4 of 9,) states that 

deg {Xi- ...■Xr)^Y[n,. 

1=1 

A refined version of Bezout's theorem (Proposition 9.1.2 of [S]) ties these results 
together with an extended Bezout formula. 

Proposition 2.4. Suppose Xi D. . . HXr consists of a connected component Z and 
a finite set S. For p G S let rUp = i{p,Xi ■ . . . ■ Xr',P^) denote the intersection 
multiplicity of p in Xi ■ . . . ■ Xr- Then deg {Xi ■ . . . ■ Xr)^ + X]pe5 "^p ~ 01=1 

If ni, . . . ,nr and X)pgs "^p known, then we can use Proposition l2.4l to solve 
for deg {Xi ■ . . . ■ X^)^ . If Z is smooth. Corollary 12.31 provides a linear relation 
among the degrees of the Chern classes of Z. We summarize these observations in 
the following theorem implicit in Fulton [S] : 
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Theorem 2.5. Let Fi,...,Fr be homogeneous forms corresponding to hypersur- 
faces Xi, . . . ,Xr C P''. Let ni = degXi, let ak he the k*'^ elementary symmetric 
polynomial in ni, . . . , n,. and let A = [oq, . . . a„] with = "^^^q {~^y {^~^^)<^7i-k-i- 
Let Z be a smooth connected n-dimensional scheme with Chern classes Cq, . • • ,c„ 
and let C = [degco, • • • ,degc„]. // the subscheme of ¥^ defined by the ideal 
{Fi,...,Fr) is a disjoint union of Z and a (possibly empty) zero-scheme S, then 
A ■ C = ar - deg S. 

In order to compute the degrees of the Chern classes of an n-dimensional smooth 
variety, it is enough to determine n+1 independent linear relations that they satisfy. 
Theorem 12.51 provides a mechanism for producing the linear relations provided we 
are able to find a sufficient number of r-tuples of homogeneous forms that satisfy 
the conditions of Theorem l2 . 51 and provided we have a method for computing deg S. 
In the next two sections, we present a collection of tools that allow us to complete 
the latter task. 

2.2. Homotopy Continuation. In homotopy continuation, a polynomial ideal, /, 
is cast as a member of a parameterized family of polynomial ideals one of which has 
known isolated solutions. Each of the known isolated solutions is tracked through a 
predictor/corrector method to a point which lies numerically close to the algebraic 
set V{I) determined by /. These points can then be refined to lie within a prescribed 
tolerance of V{L). An introduction to general continuation methods can be found 
in [T]. Through the basic algorithms of numerical algebraic geometry [35], from / 
it is then possible to produce a collection of subsets of points such that: 

• The subsets are in one to one correspondence with the irreducible compo- 
nents of the algebraic set V{I). 

• The points in a subset all lie within a prescribed tolerance of the irreducible 
component to which it corresponds. 

• The number of points in the subset is the same as the degree of the irre- 
ducible component. 

• The subset is a numerical approximation of the intersection of the irre- 
ducible component with a known linear space of complementary dimension. 

Note that the decomposition above allows one to determine both the dimension 
and degree of each algebraic variety appearing in the decomposition of an algebraic 
set. 

2.3. Degree of a Zero Dimensional Scheme. To each irreducible component of 
an algebraic set, one can use the defining set of polynomials to attach a positive inte- 
ger, called the multiplicity, that determines roughly how many times the component 
should be counted in a computation. In [6] , Dayton and Zeng study the multiplicity 
inspired, to a large degree, by Macaulay's inverse systems approach. They provide 
an algorithm which yields as output the multiplicity of isolated solutions. This is 
essentially done by counting how many partial derivatives of the polynomials are 
forced to be zero. In [B], an alternate approach is presented, inspired, to a large 
degree, by certain Grobner basis calculations coupled with a fundamental result of 
Bayer and Stillman on regularity [4]. Both of these algorithms have been imple- 
mented in Bertini Thus, there is an implemented algorithm in place that allows 
one to determine the degree of a zero-dimensional scheme. 
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3. Algorithms 

In this section, we present the pseudocode for three algorithms. The first 
algorithm computes the equivalence of a connected scheme Z in an intersection 
product. The second algorithm determines a linear relation satisfied by the degrees 
of the Chern classes of a smooth connected scheme. The third algorithm computes 
the degrees of the Chern classes of a smooth connected scheme. For each of the 
algorithms, the input is a set of non-zero homogeneous generators for an ideal 
7. It is an assumption of the algorithms that the scheme determined by / is a 
disjoint union of Z and a zero-scheme S. The second and third algorithms have 
the additional assumption that Z is smooth. 

Algorithm 1. Equivalence-of-Z {{Fi,F2, . . . , F^}; D) 

Input: A set of r homogeneous polynomials {Fi, F2, . . . , Fr} C C[zo, zi, . . . , Zr]- 
The polynomials should generate an ideal I whose corresponding scheme is the 
disjoint union of a connected scheme Z and a possibly empty zero-scheme S. 

Output: D = deg(Xi • ... ■ Xj.)^ where Xi is the hypersurface corresponding to 

Algorithm: 

Determine the support of S. 

Determine the multiplicity of each point in the support of S. 

Add up the multiplicities of the points in the support of S and store in fig. 

For each i, determine the degree of F^ and store in rij. 

Compute T = Y[i n-i- 

Compute T — fig and store in D. 

Algorithm 2. Linear-Relatiori-On-Chern-Numbers ({Fi, F2, . . . , Fr}; A, D) 

Input: A set of r homogeneous polynomials {Fi, F2, . . . Fr} C C[zo, zi, . . . , Zr]. 
The polynomials should generate an ideal I whose corresponding scheme is the 
disjoint union of Z and S, where Z is a smooth connected scheme and 5 is a 
possibly empty zero-scheme. 

Output: A = [ao, . . . ,a„] and D £ Z, where n denotes the dimension of Z. If 
Cj denotes the i*'' Chern class of Z and C = [degco, . . . , degc„] then the linear 
relation is A • C = £>. 

Algorithm: 

Determine the dimension of Z and store in n. 
Compute Equivalence-of-Z ({Fi, F2, . . . , Fr}; D). 

Compute the elementary symmetric functions ctq, (Ti, . . . , cr„ of ni, . . . , n^. 
Compute A = [ao, ai, . . . , a„] where ak = Er=o'(~l)*('^IV"-fc-»- 
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Algorithm 3. Chern_Numhers ({Fi, F2, . . . , Ft}] C) 

Input: A set of t homogeneous polynomials {i^i, ^2, • ■ ■ -Ft} C C[2;o, zi, . . . , ^r]- 
The polynomials should generate an ideal / whose corresponding scheme is the 
disjoint union of Z and S*, where Z is a smooth connected scheme and 5 is a 
possibly empty zero-scheme. 

Output: C = [degco, . . . , degc„] where n denotes the dimension of Z and Ci 
denotes the i*'' Chern class of Z. 

Algorithm: 

Determine the dimension of Z and store in n. 

Determine the degrees of i^i , . . . , Ff , and store the maximal degree as b. 
Set Ic — b for c=l,...,7i+l. 
For i = 1 to n + 1 

Choose r random elements from / of degrees li,l2, ■ ■ ■ Jr- 

Store the random elements as Gi, G2, ■ • ■ , G^. 

Compute Linear_Relation_On_Chern-Numbers ({Gi, G2, . . . , Gr}; Ai, Di). 
Let li = li + 1. 
Next i 

Build the (n + 1) x (n + 1) matrix M whose i*^ row is A^. 
Build D = [i?i,...,i?„+i]. 

Set up the linear system MC = D, with C and D as column vectors. 
Solve the linear system for C. 

There are two potential problems that should be addressed in the use of Algo- 
rithm 3. The first potential problem is whether the matrix M has full rank. The 
second potential problem is whether the random r-tuples Gi, G2, . . . , Gr, produced 
in such a simple manner from /, satisfy the input requirements of Algorithm 2. 
These two questions are answered in Proposition 13.11 and Corollary 13. 4[ respec- 
tively. 

Proposition 3.1. The matrix M defined in Algorithm 3 satisfies det(M) = ±1. 

Proof. View M as a function of the maximal degree b. Let Ui be the i*'' elementary 
symmetric polynomial in r variables bi, . . . ,br and let 

a{ = ai{bi + 1, 62 + 1, . . . , 6j 4- 1, . . . , br). 

Define r/ G by r/ = (t^(6, 6, . . . , 6), rO = a,(6, b,...,b). Define a + x + 
matrix N{b) by 

rO 1 \ 
rl 1 



N = 



\ _n n 

\ ~n 'n-l 



Tf 1 J 
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Observe first that M can be factored as M = NR, where 



(I ... \ 

r + 1 ... 



R 



-(r + l) 1 ... 



= 0. 



Since det(i?) = 1, we have det(M) = det(7V). 

One checks that det(A^(0)) = ±1. We shall show that the determinant of N{b) 
does not depend on b by showing that 

d(det(iV)) 
db 

We shall use the following fact: 

Let Ci be the i*'* column of A^, N — (Cq, Ci, . . . , C„). Then we have that 

djdetjN)) ^ ^ dC 

— 2^ aet(_Oo, Oi, . . . , Oj_i, Oj+i, . . . , 0„J, 

but each term in the above sum is zero since = (r— n+z+l)Ci+i for < i < n—1 
and = 0. □ 

Let / be an ideal defining a scheme X <zW and let Ix be the homogeneous ideal 
of X. For d G N, we say that X is cut out scheme theoretically in degree d if the 
saturation of the ideal generated by I{d) is equal to Ix- 

Proposition 3.2. Let I be an ideal defining a scheme in P'' that is the disjoint 
union of a smooth connected scheme Z and a possibly empty zero-scheme S . Let 
Gi, . . . , Gfc G /, 1 < fc < r, generate an ideal that defines a scheme X , such that 
the singular locus Xs\ng satisfies codim(Xsing \ S) > k, and codim(X \ Z) > k. 
Let d N be such that Z U S is cut out scheme theoretically in degree d and fix 
integers nk+i, . . . ,nr with Ui > d for all i. If Gk+i, . . . ,Gr are general forms in I 
with degGi ~ Ui, then the ideal J = (Gi, G2, . . . , Gr) defines the disjoint union of 
Z and a possibly empty zero-scheme S' . If S is non-singular or empty, then S' is 
non-singular. 

Proof. Suppose k ^ r. A generic form Gk+i € / of degree Uk+i cuts down the 
dimension of every component of X\(ZUS) as well as every component of Xsing\S. 
By induction, the scheme Y defined by J is a disjoint union of Z and a zero-scheme 
5", and Y is nonsingular away from S. In particular Y is non-singular on Z. □ 



Corollary 3.3. Suppose I and d are as in Proposition [37^ // Gi, . . . , Gk & I is 

a regular sequence in I such that the scheme X defined by the ideal (Gi, . . . , Gk) 
is generically smooth, then for general forms Gfc+i, . . . ,Gr with degG; = Ui > d, 
the ideal J ~ (Gi, G2, . . . , Gr) defines the disjoint union of Z and a possibly empty 



zero-scheme S' . 



Proof. Every component of X has codimension k and is of multiplicity 1. We may 
therefore apply Proposition 13.21 □ 
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The following corollary implies that for generic choices of the forms Gi , . . . , Gr 
in Algorithm 3, the assumptions on the input of Algorithm 2 are satisfied. 

Corollary 3.4. Let I be as in Proposition \3.S\ and let {Fi, . . . ,Ff} he a set 
of non-zero generators of I . Put d = max{degi^i, . . . , degi<i} and fix integers 
ni, . . . with Ui > d for all i. If Gi, . . . , Gr are general forms in I with deg Gi — 
Ui, then the ideal J = (Gi, G2, • . • , Gr) defines the disjoint union of Z and a possibly 
empty zero-scheme S' . 

Proof. Observe that the ideal generated by I{d) is equal to ®e>dl{e). Since the 
ideal ®e>dl{&) and / have the same saturation, namely the homogeneous ideal of 
Z U S, the scheme Z U S is cut out scheme theoretically in degree d. □ 

For reasons of efhciency it is desirable to keep the degrees of the forms Gi , . . . , Gr 
in Algorithm 3 as low as possible. We give two examples that have bearing on the 
question of how much Corollary [33] can be strengthened. 

Example 3.5. Consider a conic Z C P^. The conic is cut out by a plane P and a 
quadric Q. The corresponding forms of degree 1 and 2 generate the homogeneous 
ideal / of Z. In this case we cannot choose the generic forms Gi,G2,G3 G / 
of degrees (3,1,1). A cubic and two planes, all three containing the conic, will 
intersect in the conic union a line, since the planes are necessarily both equal to P. 
The residue in this case is thus a line and not a finite scheme as required. 

Example 3.6. This example is along the same lines as Example 13.51 Let A be a 

2x3 matrix of general linear forms in C[xo, . . . ,X4]. Let F2, F^, F4 be the 2x2 
minors of A and let Fi be a general form of degree 3. The ideal / = (Fi, ^2,^3, F4) 
is the homogeneous ideal of a curve Z in and the minors ^2,^3, F4 define a 
surface. We thus have a minimal generating set of / of degrees (3, 2, 2, 2) but 4 
general forms of degrees (4, 2, 2, 2) will define a union of the curve Z and another 
curve, violating the input requirements of Algorithm 2. 

4. Examples 

In this section we present several examples of computations of Chern numbers 
which illustrate the algorithms of the previous section. 

4.1. Curves. We consider the case of a smooth connected curve in P*" ([S] Example 
9.1.1). If Z is a smooth curve of genus g in P*" then the first Chern class of Z is —Kz 
and deg {—Kz) — 2 ~2g. If (Fi, F2, . . . , Fr) C C[zo, zi, . . . , Zr] is a homogeneous 
ideal I whose corresponding scheme is the disjoint union of the curve Z and a 
zero-scheme 5* then CoroUar v 1 2 . 31 leads to 

deg (Xi • . . . ■ Xr)^ = (ni + . . . + nr - (r + 1)) deg Z + 2 - 2g 

where Xi denotes the hypersurface corresponding to Fi and Ui — degX^. From 
Theorem 12. 5[ we can also write the equation as 

(2) - deg5 = (ni + . . . + rv - (r + 1)) degZ + 2 - 2g. 

i 

Example 4.1. The homogeneous ideal of the twisted cubic curve in P^ is / = 
{x"^ — wy^y"^ — xz,wz — xy) C C[w,x,y,2]. It is well known that this curve has degree 
3 and genus 0. If we choose Fi,F2,F3 G / of degrees (2,2,2) then the numerical 
irreducible decomposition implemented in [5] determines that the corresponding 
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scheme consists of a degree 3 curve and no additional points {S is empty) . Equation 
([21) gives the relation 

2-2-2-0 = (2 + 2 + 2- (3 + 1)) - 3 + 2-25 

which we can solve to get g = 0. 

If we choose Fi, F2, F3 e / of degrees (2, 2, 3) then we obtain a degree 3 curve 
and one additional point of multiplicity 1. Equation ^ gives the relation 

2-2-3-l = (2 + 2 + 3-(3 + l))-3 + 2-25 

again leading to g = 0. 

If we did not know the degree of Z, the two computations in this example would 
yield 

2 - 2 • 2 - = (2 + 2 + 2- (3 + l))dcgZ + 2- 2g, 

2 - 2 • 3 - 1 = (2 + 2 + 3- (3 + l))dcgZ + 2- 2g. 

The unique solution to these two equations is deg Z = 3 and 5 = 0. 

4.2. Surfaces. Now consider the case where Z is a smooth connected surface in P'' 
([9] Example 9.1.5). Let / = (Fi, _F2, • ■ • , Fr) C C[zo, zi, . . . , Zr] be a homogeneous 
ideal. If the scheme corresponding to / is the disjoint union of the surface Z and a 
zero-scheme S then CoroUarv 12.31 and Theorem 1 2 . 51 lead to the equation 

(3) Y[ ""-i ~ "^"^g S = deg {Xi ■ .. . ■ Xr)^ = 02 deg C2 + oi deg ci + uq deg Z, 

i 

where 

r r 

02 = 1, ai — Hi — (r + 1) and qq — Uifij — (r + 1) Ui + 

i—l i<j i—1 

Example 4.2. In characteristic zero and up to standard modifications, the 
Horrocks-Mumford bundle E is the only known indecomposable rank 2 vector bun- 
dle on P'* [16]. For a general section s of E, the zero set V{s) is a smooth surface in 
P^ called a Horrocks-Mumford surface. We shall compute the degree of a smooth 
Horrocks-Mumford surface Z as well as the degrees of its Chern classes ci and C2- 
It is known that Z is an Abelian surface of degree 10 [TB]. Hence degco = 10 
and ci = C2 = 0. The homogeneous ideal J of is generated by three quintics 
and fifteen sextics and Z is cut out scheme theoretically by three quintics and one 
sextic. Generators for / can be found by a variety of methods, e.g. by a Beilinson 
monad [7] . In this example we do not follow Algorithm 3 in detail but the example 
is covered by Proposition 13. 21 

The following table shows the number of solutions resulting from zero- 
dimensional, 128-bit precision runs in [2] on 4 random elements of / for various 
choices of degrees. We let (ni, 712, «3, "4) denote the degrees of these elements. 



(rii,n2,n3,n4) 


non-singular points 


singular points 


(5,5,5,6) 





750 


(5,5,6,6) 


40 


860 


(5,6,6,6) 


100 


980 



In each of these computations the zero set consists of the Abelian surface together 
with a finite set S of points of multiplicity one. In each row of the table the entry 
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in the middle column is the cardinality of S and the entry in the rightmost column 
is the Bezout number Y[t=i ^« minus the entry in the middle column, i.e. 

nin2n3?M — rUp. 

pes 

By Proposition 12. 4i this is also the degree of the equivalence of Z. Evaluating Q 
in each case gives: 

deg C2 + 16 deg ci + 75 deg Z = 750 
deg C2 + 17 deg ci + 86 deg Z = 860 
deg C2 + 18 deg ci + 98 deg Z = 980. 

This system has the unique solution degZ = 10, degci = degC2 = 0. 

4.3. Higher Dimensional Varieties. If Z is a smooth n-dimensional variety in 
P'' then the algorithm proceeds in a similar manner to the curve and surface cases. 
We use homotopy continuation to determine the equivalence of Z in various inter- 
section products. This combines with the formulas appearing in CoroUarv 12.31 and 
Theorem 12.51 to produce a linear system involving the degrees of the Chern classes 
of Z. Finally, we solve the linear system to determine these Chern numbers. 

Example 4.3. Let / be the ideal defined by the 4x4 minors of a 4 x 5 matrix of 
general linear forms in C[xq, xi, . . . ^x^] and let Z be the corresponding threefold 
in P^. The following table shows the number of solutions resulting from zero- 
dimensional, 128-bit precision runs in [2] on 5 random elements of / for various 
choices of degrees. We let (7^1,712,713,714,715) denote the degrees of these elements. 



(7^l,?^2,?^3,'^4,?^5) 


non-singular points 


singular points 


(4,4,4,4,4) 





1024 


(4,4,4,4,5) 


1 


1279 


(4,4,4,5,5) 


6 


1594 


(4,4,5,5,5) 


21 


1979 



In each of these computations the zero set consists of the threefold together with 
a finite set S of points of multiplicity one. The formulas from CoroUarv 12.31 and 
Theorem 12.51 lead to: 

deg C3 + 14 deg C2 -t- 61 deg ci + 44 deg Z = 1024 
deg C3 + 15 deg C2 -t- 71 deg a + 65 deg Z = 1279 
deg C3 + 16 deg C2 + 82 deg ci + 92 deg Z = 1594 
degC3 + 17degC2 -H94degci + 126 deg Z = 1979. 

The unique solution is degZ = 10, degci = 0, degC2 = 45, degcs = —46. 

Example 4.4. Let V be the image of the Segre embedding i : P^ x P^ ^ P*". Let 
Z be the intersection of V with a general hyperplane and let / be the homogeneous 
ideal of Z . The following table shows the number of solutions resulting from zero- 
dimensional, 128-bit precision runs in 2j on 6 random elements of / for various 
choices of the degrees (tii, . . . , ng) of these elements. 



(711,712,713,714,715,716) 


non-singular points 


singular points 


(2,2,2,2,2,2) 





64 


(2,2,2,2,2,3) 





96 


(2,2,2,2,3,3) 


2 


142 


(2,2,2,3,3,3) 


10 


206 
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In each of these computations the zero set of the 6 random elements consists of the 
threefold together with a finite set of points of multiplicity one. The formulas from 
Corollary 12.31 and Theorem 1 2 . 51 lead to: 

deg C3 + 5 deg C2 + 4 deg ci — 8 deg Z = 64 
deg C3 + 6 deg C2 + 7 deg ci — 10 deg Z — 96 
degC3 + 7degC2 + lldegci - lldegZ = 142 
deg C3 + 8 deg C2 + 16 deg ci - 10 deg Z = 206. 

The unique solution is dcgZ = 4, degci = 10, degC2 = 10, degca = 6. 

5. Conclusion 

The results of this paper demonstrate the viability of computing Chern numbers 
of smooth varieties through numerical homotopy continuation. Homotopy continu- 
ation via square systems is a natural venue through which to compute the equiva- 
lence of a scheme in the context of an important generalized Bezout's theorem from 
intersection theory. It should be noted that the algorithms of this paper could be 
implemented in a purely symbolic setting as well. The advantages of the numerical 
approach is that the algorithms work equally well whether the generators are sparse 
or dense and whether they have rational, algebraic or transcendental coefficients. 
In addition, meaning can often be attached to the computations in situations where 
the generators have inaccuracies in their coefficients. Finally, homotopy continua- 
tion algorithms parallelize well allowing such algorithms to take full advantage of 
multi-processor machines. These features suggest the complementary nature of the 
approach to purely symbolic methods. 

In future work, the authors intend to extend the approach to take further ad- 
vantage of the ideas of intersection theory. In particular, intersection theory for 
non-square systems leads to a method for computing intersection numbers of Chern 
classes. 
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